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Abstract 

Between 2004 and 2007 the instruments of the CASSINI spacecraft, orbiting 
within the Saturn system, discovered dark patches in the polar regions of 
Titan. These features are interpreted as hydrocarbon lakes and seas with 
ethane and methane identified as the main compounds. In this context, we 
have developed a lake- atmosphere equilibrium model allowing the determina- 
tion of the chemical composition of these liquid areas present on Titan. The 
model is based on uncertain thermodynamic data and precipitation rates of 
organic species predicted to be present in the lakes and seas that are subject 
to spatial and temporal variations. Here we explore and discuss the influence 
of these uncertainties and variations. The errors and uncertainties relevant 
to thermodynamic data are simulated via Monte-Carlo simulations. Global 
Circulation Models (GCM) are also employed in order to investigate the pos- 
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sibility of chemical asymmetry between the south and the north poles, due to 
differences in precipitation rates. We find that mole fractions of compounds 
in the liquid phase have a high sensitivity to thermodynamic data used as 
inputs, in particular molar volumes and enthalpies of vaporization. When 
we combine all considered uncertainties, the ranges of obtained mole frac- 
tions are rather large (up to ~ 8500%) but the distributions of values are 
narrow. The relative standard deviations remain between 10% and ~ 300% 
depending on the compound considered. Compared to other sources of un- 
certainties and variability, deviation caused by surface pressure variations 
are clearly negligible, remaining of the order of a few percent up to ~ 20%. 
Moreover no significant difference is found between the composition of lakes 
located in north and south poles. Because the theory of regular solutions em- 
ployed here is sensitive to thermodynamic data and is not suitable for polar 
molecules such as HCN and CH3CN, our work strongly underlines the need 
for experimental simulations and the improvement of Titan's atmospheric 
models. 

Keywords: planets and satellites: individual: Titan - planets and 
satellites: general - solar system: general 



1. Introduction 



The surface of Saturn's haze-shrouded moon Titan had long been pro- 
posed to be at least partly hidden by oceans or s eas, on the basis of the sta- 



bility of liquid methane an d ethane at the surface (jFlasarl . ll983l : lLunine et al. 



19831 : iLorenz et al.l . |2003| ). The presence of a global ocean on Titan was ex- 



clude d from ground-based radar observations in the mid 1990s (IMuhleman et al. 



I995I ). In mid 2006, dark, lake-like features of a rang e of sizes were detec ted at 
Titan's north polar region by the Cassini RADAR (jStofan et al.l . 120071 ). The 
chemical composition of these lakes remains, however, poorly determined. 
Spectra of the southern hemisphere lake Ontario Lacus have been obtained 
by the Visual and Infrared Mapping Spectrometer (VIMS) a board Cassini 
but t he only species that has been firmly identified is C2H6 ( IBrown et al. 
20081 ). The difficulty in determining the composition of the lakes is essen- 



tially due to the presence of a large atmospheric fraction of CH4 that impedes 
this molecule's identification in the liquid phase present on the surface, ir- 
respective of the value of its mole fraction. However, methane is indirectly 
inferred in Ontario Lacus by the secular decline of the lake extent over the 
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Titan summer ( iHayes et al.l . |2010| ) and the observation of troposphere clouds, 
which must be methane, coincident with sur face darkening over the southern 
pole during the summer (ITurtle et al.l . 120071 ). Because the detection of other 
compounds in the lakes of Titan remains challenging in the absence of in 
situ measurements, the only way to get a good estimate of their chemical 
composition is to develop and utilize a thermodynamic model based on the- 
oretical calculations and laboratory data. Several models investigating the 
influence of photochemistry and the atmospheric composition on the chemi- 
cal composition of putative hydrocarbon oceans or seas f ormed on the surface 
of Titan, have been elaborated in the pre- Cassini years (ILunine et al.l . Il983 : 
Dubouloz et al.l . Il989l : iMcKay et al.l . Il993l : iTokanol . 120051 ) . These models sug- 
gested that the liquid phase existing on Titan contains a mixture made from 
C2H6, CH4 and N2, and a large number of dissolved minor species. 



On the other hand, the Cassini- Huygens measurements have improved our 
knowledge of the structure and composition of Titan's atmosphere. In partic- 
ular, the Gas Chromatograph Mass Spectrometer (GCMS) aboard Huygens 
and the Cassini Composite Infrar ed Spectrometer (CIR S) provided new at- 
mospheric mole fraction data (see iNiemann et al.l . l2005l . and Table 1). More- 
over, near-surface brightness temperatures and corresponding estimates for 
physical temperatures in the high latitudes a t which numerous lakes are found 
have now been determined (jjennings et al.l . |2009|). These atmospheric and 
surface conditions have been recently used to re compute the solubili ties of the 



different compounds in the hydrocarbon lakes (ICordier et al.l . l2009l hereafter 
C09). The same model has also been employed to explore the possibility of 
noble gas trapping in the lakes of Titan in order to proy i de an attempt of 



explanation of their atmospheric depletion (ICordier et al.l . I2OIOI ) 



The assumptions considered by C09 are similar to those made by lDubouloz et al. 



( 1l989l ) (hereafter DUB89): in both cases, lakes are considered as nonideal 
solutions in thermodynamic equilibrium with the atmosphere. However, nei- 
ther DUB89 nor C09 have taken into consideration the influence of uncer- 
tainties on the data used as inputs in their models. Indeed, some thermo- 
dynamic data are measured at much higher temperature and extrapolated 
down to temperatures relevant to Titan's conditions. Precipitation rates are 
also supposed to vary with respect to latitude, longitude and time. In this 
work, we investigate the influence of thermodynamic uncertainties, and in a 
lesser extent, the geographic influence of the variation of precipitations on 
the lakes composition. In the latter case, we restrict our study to a supposed 
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north/south poles asymmetry in chemical composition. 

In Section H] we detail our lake-atmosphere equilibrium model. Section [3] 
is dedicated to the study of the influence of uncertainties on thermodynamic 
data (vapor pressures, molar volumes, enthalpies of vaporization and parame- 
ters of interaction) on the resulting lakes composition. In Sect. HI simulations 
are conducted with the use of precipitation rates derived fro m a version of 
the IPSlE 2-dimensional climate model of Titan's atmosphere ( ICrespin et al 



20081 ) and allow comparison between chemical composition of south pole and 



north pole lakes. Section |5] is devoted to discussion and conclusions. 

2. Description of the lake-atmosphere equilibrium model 

Our model is based on regular solution theory and thermodynamic equi- 
librium is assumed between the liquid and the atmosphere. This equilibrium, 
which is expressed by the equality of chemical potentials, can be written as 
follows (Eq. 1 of DUB89): 

^fc -P = Tfc Pvp,k, (1) 

where P is the total pressure at Titan's surface, and respectively the 
mole fractions of the k compound in the atmosphere and in the liquid, and 
Pvp,k its vapor pressure. The activity coeflicient Tk (dimensionl ess) of the k 



compo und is given by (frame of the regular solution theory - see iPoling et al 



mom -. 

RTlnT, = Vrn,k J2 SZ^^^'^ - ^^^^ (2) 



« 3 

where 



Aj = {Si - 5jY + 2/ij 5i5j (3) 

and 

$^=XiV;„,^/5Z^iKnJ. (4) 



Institut Pierre-Simon Laplace 
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5i ((J.m ^)^/^) is the Hildebrand's solubility parameter of the ith compound. 
The value of this parameter is given by: 



, AHjj i — RT . , 



where AH^i (J mol~^) is the enthalpy of vaporization and Vm^i (m^ mol~^) the 
molar volume. A Si represents a measure of the molecular cohesion energy 
of the pure component i. It depends on the nature and the strength of 
intermolecular forces (hydrogen bond, ...) between molecules of the same 
species. In general, two components i and j with 6i and Sj presenting close 
values, have a high solubility. Beside this, the /jj's parameters represent the 
effects of interactions between molecules of different species. These Zj^'s are 
empirically determined and are generally poorly known. The situation Vi, j: 
Si = Sj and lij = corresponds to all activity coefficient equal to one, in 
other words this is an ideal solution in which all intermolecular forces are 
negligible. 

Our model also allows us to estimate the mole fraction of each solid 
precipitate that is dissolved in the lakes of Titan. To this end, we calculate 
the saturation mole fraction|§ Xj^^at of the compound i, which is given by (Eq. 
7ofDUB89): 

\n{r,Xi,sat) = {AHjRTJil - TJT), (6) 

where is the component's melting temperature and AH^ its enthalpy of 
fusion. Our calculation procedure is then as follows: 

1. The unknown Xj's and Fj's are computed via the Newton-Raphson 
method. 

2. Once the Xj's have been determined, the Xi^sat^ are in turn calculated 
and compared to the Xj's for each species. If for i compound we get 
Xi,sat < Xi, then we fix Xi = Xi^^at- 

3. We get new values of Xj's and Xj^sat's via the resolution of the nonlinear 
system. 



^The saturation mole fraction of the compound i corresponds to the maximum mole 
fraction of t in the liquid form. Above this value, the i material in excess remains in solid 
form. 
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4. The iterations are continued until we get a difference between Xi^sat 
and Xi lower than 10~^, value for which the numerical inaccuracy is 
clearly negligible compared to other sources of uncertainties. 

The known Fj's are given in Table [TJ The precipitation rates r^'s represent 
the number of molecules of a given species, reaching the surface of Titan by 
unit of time and by unit of surface (molecules m~^ s~^). The r, ;' s used in 



C09 were derived f rom t he photochemical models of iLavvas et al.l (j2008aU bl) 



and IVuitton et al.l (120081 ) . and correspond to the main products of CH4 and 
N2 photolysis. These rates allow us to express each i compound that falls 
from the atmosphere in the form 

Xi = X Xc,H, (7) 

We also ensure that YliXi = 1 and = 1. In this way, we get 15 un- 

knowns and 15 equations, allowing the system to be solved. The thermod 
namic data used in our calculations derive mainly from the NIST databasi 
As discussed in the following section, these data are often not well known and 
the large uncertainties associated to their determination may induce strong 
variations of the lakes chemical composition. 

3. Uncertainties due to thermodynamic data 

In the calculations of C09, the thermodynamic data (vapor pressures, 
molar volumes, enthalpies of vaporization and /ij's) have been set to their 
nominal values in the lake-atmosphere equilibrium model. However, each of 
these nominal values is accompanied by a given "deviation" or "error" and 
the consideration of the full range of possibilities for these thermodynamic 
data may strongly alter the lakes composition compared to the one calcu- 
lated by C09. In order to investigate up to which point the composition of 
these lakes may depart from the one of C09, we use here a M onte-Carlo nu- 



merical method (see for instance [Metropolis and Ulaml . Il949[ ). allowing us to 



perform error simulations for vapor pressures, molar volumes, enthalpies of 
vaporization and parameters of interaction Zjj's. A first set of computations 
consists in calculating the composition of lakes from numbers randomly cho- 
sen within the range of possible values attributed to a set of thermodynamic 
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data. For each mole fraction Xj, the minimum Xi^min, the maximum Xi^^nax 
and the average value Xi are recorded. The procedure is repeated 10,000 
times for each set of thermodynamic data. The choice of the total number of 
Monte-Carlo iterations is a compromise that has been fixed to get a statisti- 
cally significant population while maintaining reasonable computation time. 
Simulations with 5,000 and 20,000 iterations do not give significantly differ- 
ent results. Additionally we perform a simulation addressing the case where 
all sets of thermodynamic data are simultaneously considered with synthetic 
errors. This procedure is also repeated 10,000 times. 

Note that the enthalpies of melting AH^ are not considered in our inves- 
tigation because i) they seem to be reasonably well known compared to other 
thermodynamic quantities and ii) they only play a role in case of saturation, 
e.g here with HCN. For this compound the measurement, p rovided by the 
NIST database, comes from iGiauque and RuehrweinI ( 119391 ) with an accu- 
racy of about 10~^. Such a level of uncertainty is clearly negligible compared 
to other sources, allowing us to keep these thermodynamic quantities out of 
our study. 

3.1. Influence of vapor pressure uncertainties 

Vapor pressures of species, for which Eq. [1] is written, are taken from the 
NIST database in the form of an Antoine's law in the cases of N2, CH4, Ar 



and C2H6 or from a vapor pressure law given by iLidd ( 1l974l ) in the case of 
CO0. In general the domains of validity of Antoine's laws used in this work 
include the range of temperatures relevant for Titan's lakes [i.e. 90 ± 3 K). 
For instance, in the cases of CH4 and C2H6, the lower boundaries are 90.99 K 
and 91.33 K respectively, implying moderate extrapolations for temperatures 
slightly below ~ 90 K. 

On the other hand, evaluating the accuracy of Antoine's equations brought 
by the NIST database is not straightforward. To do so, we have first con- 
sidered the c ase of N2 for which NIST maintainers derived an Antoine's 



equation from lEdejer and ThodosI (119671 ). These authors published a Frost 



Kalkwarf equation based on 180 experimental vapor pre ssure measurements 



derived from 13 references. The equation obtained by lEdejer and Thodos 



(119671 ) reproduces the experimental measurements with a deviation rang- 



ing between 0.13% and 2.04%. Comparing the vapor pressure computed 



^CO data are unavailable in the NIST. 
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with t he NIST Antoine's equation with the one given by Edejer and Thodoi 
(Il967l ). we found differences reaching ~ 10% for the lowest temperatures {i.e 
around 67 K) and ~ 1% for temperatures close to 90 K. Following a similar 
approach, we compared pr essure computed with the NIST Antoine's equa- 
tion and original data from lCarruth and Kobayashil (jl973[ ). We also ma de a 
comparison between our own fit and pressure data given by iLidd (1X9741 ) for 
carbon monoxide. For the relevant temperature domain, the deviation for 
C2H6 remains between 0.1% and 1%, while the vapor pressure of CO reaches 
a difference of about 9%. Consequently we have fixed the maximum errors on 
vapor pressures, for all relevant species (i.e. N2, CH4, Ar, CO and C2H6) to 
±10% relative to previously used values (see C09). This range should bracket 
all the vapor pressures expected for each compound. This approach allows 
us to explore a wide range of possibilities, including combinations which do 
not correspond to physical reality. In this sense, results corresponding to 
extreme deviations should be regarded as unlikely cases. 

Table [2] gathers our results which are quantified by A p^^^ = (X^ax ~ 
Xmm)/X and the relative standard deviation a* (both expressed in percent- 
age). Ap^^pmeasures the total spread (over 10,000 computations of chemical 
compositions) of mole fraction values for a given species, including the results 
of the most unlikely combinations of synthetic errors. The relative standard 

deviation a*{i) = ^Jxf - X^jXi (the upper bar denotes the average value 
over 10,000 computations) shows how much variation or "dispersion" there is 
from the "average" Xi. Table|2]shows that a*{i)^s differ strongly from Ap^^p's. 
This feature corresponds to the signature of extremely narrow distributions 
of values around the average ones and can be explained by the non-linearity 
of the equations of our model. Indeed, synthetic errors are chosen with a 
uniform distribution but the resulting distribution of mole fractions is het- 
erogeneous. The shapes of Xj's distributions are shown in subsection 13.51 in 
which errors for all thermodynamic inputs are taken into consideration. 

Errors on vapor pressure mainly affect mole fractions of species for which 
Equation [1] is written. This behavior is not surprising as Equation [1] contains 
explicitly the vapor pressure. The entire set of equations being coupled, 
even a variation of one vapor pressure affects the mole fractions of all the 
other species determined with our model. Note that the case of ethane is 
particular because its atmospheric mole fraction yc2H6 is an unknown of our 
mathematical problem, while atmospheric abundances of N2, CH4, Ar and 
CO are fixed by the observations. Compared to other compounds belonging 
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to precipitated species HCN shows a relatively high Ap^^ (c"*)- In this case, 
this is also due to the use of an equation (Equation E]) where vapor pressures 
play a role via Fj. Average values Xi differ slightly from previous results 
(C09). 

3.2. Influence of molar volume uncertainties 



Molar volumes have been estimated v ia Rackett's metho d (see lPoling et al. 



20071 ). Table 4-11 and pages 4.36-4.37 of IPohng et aP fl2007h present compar 



isons between measured molar volumes and estimated ones for some organic 
compounds at the boiling temperature. These comparisons show that maxi- 
mum deviations typically reach the levels of a few percent. We then adopted 
a maximum "error" of ±10% for Monte- Carlo simulations only applied to 
molar volumes. These simulations allow a sensitivity comparison with those 
performed for vapor pressure. Resulting Ay^and cr^^are displayed in Table |2l 
Similarly to simulations related to vapor pressure, large differences be- 
tween Ay^and Uy^indicate a very narrow spread of mole fractions in lakes. 
As would be expected, species for which thermodynamic Equations [T] are 
explicitly used show the highest deviations. As also shown by Table [21 a gen- 
eral trend is that mole fractions appear to be much more sensitive to molar 
volume than to vapor pressures. 

3.3. Influence of enthalpies of vaporization uncertainties 

The NIST and the literature provide numerous interpolation formulae for 
enthalpies of vaporization. Getting reliable estimates of their actual accu- 
racy is not easy because the domains of validity of these formulae do not 
often include the ground temperature of Titan. For instance for C2H6, CaHg 
and C4H8 extrapolations over about 100 K are required. However we have 
performed error estimates only for methane, ethane and argon. In the cases 
of methane and ethane, we compared the ent halpies of vaporization give n 
by the NIST database (originally published by Maier and Svobqda ( IQSsI l) 



to those computed with the equations provided by ISomayajuld (119881 ). For 
methane, in the temperature range of interest for the surface on Titan, we 
obtained differences lower than 1%. In the case of ethane, these differences 
are muc h more impor t ant a nd lie between 26% and 30%. Considering the 



work of iTegeler et al.l f ll999l ). we estimate an internal uncertainty of about 
1% for the enthalpy of vaporization of argon. 

Again, in order to be consistent with others Monte-Carlo simulations, we 
fixed the maximum "error" on the enthalpy of vaporization of each species 
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to ±10%, a value which is well within the range of uncertainties found from 
comparisons. The results are displayed in Tableland show that the induced 
uncertainties on mole fractions are similar or higher than those obtained for 
the molar volumes. 



3.4- Influence of Uj uncertainties 

The interaction parameters /jj's represent the interaction between molecules 
of different species and are essentially determined empirically. These param- 
eters are fixed to zero in the case of interactions between the same molecules 
(Vi la = 0). In principle these 1^ depend on temperature, however for typical 
nonpolar mixtures over a mode s t ran ge of temperature, that dependence is 



usually small (see iPoling et al.l . 120071 ). As stated by DUB89, they are un- 



known in many situations. In this work, as 13 species are taken into account 
in the liquid phase, we need to know 156 parameters (13 x 13 — 13 = 156; 
the 13 la being set to zero) and testing the possible influence of each of them 
in our system does not really make sense. Giv en the fact that th e values of 



kj range between ~ 0.02 (DUB89) and 0.09 flPohng et al.l . 120071 ). we have 
performed Monte-Carlo simulations with these parameters set randomly be- 
tween and 0.10 (except la = 0). The sensitivity of mole fractions to k/s 
is presented in Table |2] and appears lower than in the cases of molar volume 
and enthalpy of vaporization. 

3. 5. Combination of all thermodynamic uncertainties 

Here we have combined all sources of uncertainties {i.e. errors on Pvap, Vm, 
AHyap and kj) and the results are represented by Aau and in Table |2l 
The combination of uncertainties on thermodynamic data can induce mole 
fraction fluctuations up to a factor of ~ 100 for Aau and slightly lower than 
~ 4 if we consider the relative standard deviation a^w- The distribution of 
mole fractions Xi is represented in the form of histograms in Figures [1] and 

El 

For each species, the range Xi^rnin ^ ^i,max has been divided into 100 
intervals. For the k-th interval, the number A^^ of mole fractions owning 
a given value has been normalized via = A'^/Apeak, where A'peak cor- 
responds to the largest N^. We stress that Apeak has a specific value for 
each compound. Distributions are very narrow and clearly asymmetric. The 
smallest abundances are limited by X = 0. HCN is a particular case for 
which the highest abundances are limited by the saturation. Curiously, the 
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HCN distribution presents a "residual tail" located at lake mole fractions 
between and ~ 0.032. 

If we consider case 1 of DUB89 {i.e. T = 92.5 K, lAr = and 
^CH4 = 0.0155), our results bracket the abundances found by these authors. 
For instance we find 7.7 x 10"^ < Xn^ < 0.039 while DUB89 got Xn^ = 0.018. 
The case of methane is similar since we find 1.8 x 10~^ < Xqu^ < 0.083 
whereas DUB89 inferred XcUi = 7.3%. This illustrates the fact that differ- 
ences between C09 and DUB89 are consistent with uncertainties caused by 
poorly known thermodynamic data. 

One could argue that the choice of a maximum deviation of ±10% is 
arbitrary, even if that level of uncertainty has been discussed in previous 
subsections. Figure [3] shows the sensitivity of a* to the adopted maximum 
error for N2, CH4, Ar and CO. As expected, the standard deviation increases 
with the value of the maximum error but this behavior appears to be non- 
linear. During a Monte-Carlo simulation, some combinations of errors yield 
to a non-convergence of the model, these occurrences corresponding more 
likely to unphysical situations and/or to an initial input (in practise the 
nominal solution for T = 90 K published in C09) in the Newton-Raphson 
algorithm which is too far from the solution of the system of the equations. In 
addition, a maximum error of ±10% in the thermodynamic data is probably 
an overestimated value even if we do not really know how these data depart 
from the "real" ones. However, we consider that the most important point 
here was to address the sensitivity of the model to the different sources of 
uncertainties. 



4. Influence of geographic variations of precipitation 

On Earth, precipitation is almost entirely water, with the rate dependent 
on location and time. In the case of Titan, the situation is more complex 
because the slow sedimentation of stratospheric aerosols to the surface is key 
to filling the lakes with the dominant photochemical byproducts of methane 
as well as less abundant species. 

Moreover, the sedimentation rates are also a function of location and time. 
In this section we restrict our study to the geographic dependence of lakes 
composition and more precisely to a possible South/No rth asymmetry as 



the d istribution of lakes seems to be itself asymmetric (see lAharonson et al 



20091). 
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As mentioned in Sect. [21 the calculations are based on p r ecipita. tion rates 
derived from a slightly improved version of iLavvas et al. teoOSaU bl) (here- 
after LAV08) models that included the atmospheric profile of Titan's at- 
mosphere measured with the Huygens Atmosphere Structure Instrument 



(iFulchignoni et al.l . |2005[ ) and some updated reaction rates. In the follow- 
ing we refer to this set of precipitation rates as LAV09. These models are 
clearly inadequate for a ge ographical study, this is why we used a 2D mod- 
els originally developed by iLebonnois et al.l ( 1200 ll ) . These authors used an 
analytic description of the meridional circulation of Titan's atmosphere to 
take advection into account in a two-dimensional photochemical model. This 
coupling between dynamics and photochemistry was subsequently improved 
by the implementation of th e same photochemic al model in the IPSL two- 



dimensional Climate Model (jCrespin et al.l . |2008| . hereafter LEB08 



Precipitates are included in our thermochemical model via the relation 
Xi = — X Xc2He- At the first sight, one could believe that multiplying a 
given Tj by an arbitrary factor yields a multiplication by of the resulting 
mole fraction X^. This is not the case because the mole fraction of ethane 
-^C2H6 depends on all the other mole fractions. Indeed Eq [1] depends on 
EqHl and both equations are solved simultaneously. Hence, the influence of 
precipitation rates on mole fractions can only be estimated with a complete 
calculation. As these rates play a role in our set of equations via the ratios 
Ti/TcaHe! we then consider these ratios instead of absolute values. For 
species available in LEBOS's models, we computed the time averaged ratios 
Tj/rcaHe; Tcsults are displayed in Fig. HI ratios deduced from LAV08 and 
LAV09 have been displayed for comparison. 

While we have noticed in LEB08 data absolute precipitation rates have 
huge latitudinal variations, the ratios Ti/rcjHe do not exhibit such steep de- 
pendence in regions located poleward of latitudes around ±60°. In equatorial 
regions, the precipitation rates computed by LEB08 can be very small, im- 
plying that the ratio rj/rcjHe probably has no great physical meaning in 
these regions. Fortunately, as hydrocarbon lakes have presumably been de- 
tected in polar regions, this observational evidence allows us to identify the 
questionable values. 

With LEB08 polar ratios Tj/rcaHe we compute two sets of mole fractions: 
one for the south pole, another for the north pole. In both cases, the temper- 
ature and the pressure have been respectively fixed to 90 K and 1.467 bar. 
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We did not find any significant differences between the composition of south 
and north lakes (i.e. differences of the order of 1%) except for CsHg (~ 30%), 
C4H10, CH3CN and CgHe (both around ~ 20%). The north pole mole frac- 
tions being systematically larger than those computed for the south pole, this 
behavior corresponds to LEB08 large ratios at north pole (see Fig. H]). These 
results have to be considered carefully because even up-to-data 2D Titan's 
atmosphere models have to be improved. For instance, LEB08 precipitation 
rates are in fact condensation rates. In this approach, when a given species 
is in an atmospheric layer where the local temperature corresponds to the 
saturation temperature of this species, then all the molecules in this layer are 
supposed to precipitate on Titan's surface. In this picture, the microphysics 
of clouds is not taken into account. Models including this microphysics for a 
lot of species have to be developed. 



5. Discussion and conclusion 



Beside vapor pressures, molar volumes and other parameters already stud- 
ied in previous sections, total pressure P and temperature T at ground level 
could have also an influence on lakes composition. 

The influence of temperature has been already discussed in C09. The vari- 
ations of ground pressure, as P appears in Eq. [H could change the thermo- 
dynamic equilibri um. Using the Cassini Synthetic Aperture Radar (SAR) 
Stiles et al.l ( l2009l ) have developed models of the topography of limited por- 
tions of the surface, finding surface heights typically in the range -1500 m 
to +1000 m, yielding a maximum altit ude difference of about 2500 m. More 
recently, radar altimetry analyzed by I Wall et al.l ( I2OIOI ) across Ontario La- 
cus and its surroundings shows a maximum amplitude of the altimetry echo 
center of mass of about ~ 500 m. This l ocal determination is co mpatible 



with IStiles et al.l ( 12009! ) work. HASI data ( iFulchignoni et al.l . |2005| ) contain 
pressure records, measured on January 14th, 2005 during the Huygens probe 
descent. Between the Huygens landing site and an altitude of 2500 m, the 
pressure ranges between 1467 hPa and 1296 hPa, i.e. a relative variation of 



10%. By means of their 3D general circulation model, iTokano and Neubauer 
( I2OO2I ) have investigated the influence of Saturn's gravitational tide on the 
atmosphere of Titan. They found that induced surface pressure variations 
remains lower than 1.5 hPa, i.e. 0.1%, which is negligible in our context. 
Table [3] summarizes our results when we change the pressure value by ±10% 
from the Huygens Huygens-derived values Puny 
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It is a valid question as to whether thermodynamic equilibrium is a rea- 
sonable assumption or not. For that purpose one can calculate the thermal 
relaxation time of a lake with a depth of H given by 



(8) 



(see the classical textbook Landau and Lifshitzl . 1987 ). with x the thermal 
diffusivity which is given by x = n/ pC-p where n (W.m~^.K~^) is the thermal 
conductivity, p (kg.m~^) the density and Cp (J.kg~^.K~'^) the mass, specific 
heat capicit y at cons t ant p ressure. We estimated n thanks to the Enskog's 



theory (see iDymondl . Il985l ). we found k ~ 0.59 W.m ^.K ^ for pure hq- 



uid ethane at T = 90 K. T h is val ue is of the same order of magnitude as 
that found by iLorenz et al.l (120101 ) (see their table 1), who got k ~ 0.25 
W.m~^.K~^. Fin ally we o btain ed r ~ 2 Titan's days for if = 1 m. As 
it can be seen in iTokand (j2005[ ) Titan's surface temperature variations in 
the polar regions are of the order of 1-2 K over a Titan's year (about 673 
Titan's dayqj), that means that -at least the first meter of the lakes- have 
enough ti me to be i n the rmal equilibrium with the atmosphere. Following the 
models of iTokanol fl2005[ ). during spring and summer the lakes are thermally 
stratified (which may imply a stratified chemical composition), but by the 
autumnal equinox convection renders the upper layer of the lakes isothermal. 



In our work we use the semi-empirical regular solution theory already ern- 



ployed in previous works ( jPubouloz et al.l . Il989l : ICordier et al.l . l2009l l2010l ) 



but this also has its limits of validity. Pa rticularly, it has been i ntrod uced 
for binary mixtures of nonpolar molecules ( Hildebrand and Scott . 19621 ) and 



Poling et al.l . 120071 chaper 



generalized to multicomponent mixtures (see 
Even though this generalization appears to be reasonable it has never been 
properly validated in a context relevant for Titan's hydrocarbon lakes. More- 
over the interaction parameters (the Uj^s) are not well known in the case of 
nonpolar molecules and are probably an inadequate formalism when polar 
components are in the solution. We recall that the mixture considered in 
this work includes two polar molecules: HCN (with a dipole moment of 2.98 
D, which could be compared to the water dipole moment of 1.85 D) and 
acetonitrile CH3CN (with a dipole moment of 3.84 D). 



^Onc Titan's day corresponds approximately to 16 terrestrial days. 
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Our work stresses the great impact on predicted composition of uncer- 
tainties in the thermodynamic inputs. In the framework of our model, this 
influence appears to be more important than abundance differences between 
north and south pole lakes, assuming a maximum "error" level of ±10% con- 
sidered in our Monte-Carlo simulations. Our computations show also that 
the influence of pressure variations is purely negligible. 

It is important to note that, in this work, we did not consider the tempo- 
ral variation of lakes chemical composition. Indeed, many phenomena could 
contribute to these variations, among which seasonal variat ions for short 



times cales and the Milankovitch cycle for longer timescales (see lAharonson et al. 
20091 ). We stress that all these phenomena involve processes of evapora- 
tion/condensation of various species (in particular CH4), which clearly rep- 
resent non-equilibrium situations. A more realistic model will have to take 
into account e nergy and mass fluxes between lakes and the atmosphere (see 
Tokand . l2005l ) and incorporate a chemical model similar to the one used in 
the present work. Future works will have to integrate these two aspects of 
modelling to provide a more accurate description of lakes evolution. 

We also underl ine the need for mor e realistic photochemistry models as 
already stated by (IHebrard et al.l . 120071 ) among others. If we concentrate on 
the properties of the hquid of the lakes themselves, two kinds of approaches 
can be considered to make progress beyond what has been done here: (1) the 
development of more accurate thermodynamic data (measured in dedicated 
experiments and/or determined by ab initio computations); (2) Titan's lakes 
in vitro simulations, in which one explicitly attempts to simulate Titan's lakes 
through a liquid hydrocarbon mixture in contact with an atmosphere in a 
laboratory chamber. The second approach is surely more relevant because 
a model is not required to apply pure thermodynamic data, but the first 
approach may be more practical in the absence of a major experimental 
effort tied to propo sed future missions to Titan like Titan Saturn System 
Mission (TSSM, se e lMatson et al.l . l2009l ) or Titan Mare Explorer (TIME, see 



Stofan et al.l . I2OIOI ). this latter being dedicated to lakes study and analysis. 
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Tabic 1: Assumed composition of Titan's atmosphere at the ground leveL 



Atmosphere 


Mole fraction 


Determination 


H2 


9.8 X 10"^ 


Huygens GCMS^"^ 


CH4 


0.0492 


Huygens GCMS^^) 


CO 


4.70 X 10-5 


Cassini CIRS^^) 


40 


4.32 X 10-5 


Huygens GCMS^'') 


N2 


0.95 


CO9W 


C2H6 


1.49 X 10-5 


CO9W 



('^^Owen & Niemann 2009; ('') Niemann et al. 2005; ("^^De Kok et al. 2007; 
('^)N2 and C2He abundances have been calculated by C09 and correspond to 
a ground temperature of 93.65 K. In Monte-Carlo simulations presented 
here, the mole fractions of N2 and C2H6 are varying (see text). 
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Table 2: Results of Monte-Carlo simulations showing the space of possible values for thermodynamic inputs explored randomly 
within a ±10% range centered on the nominal values used by COO. X„iin, X and X,„ax are recorded for each species and results 
are presented in the form of A's and relative standard deviations a* (see text). X^i^^ X and Xmax are shown only in the cases 
where the space of vapor pressures is explored. Ap^^^ (cTp^^^ ) refers to Monte-Carlo simulations for which only vapor pressures 
are affected by synthetic errors. Av'^((Ty^), Aap^^p (c^ji/^^p), and /S.i^.{a*..) refer to molar volumes, enthalpies of vaporization 
and interaction parameters Uj^ respectively. Aaii(ct^ii) are resulting uncertainties when all thermodynamic quantities are 
considered with errors. 





Y 

^min 




X 




^max 




vap 


vap 










A/.. 




Ami 


* 




(-Pvap) 




(-Pvap) 




(-Pvap) 




(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


(%) 


N2 


4.1 • 10' 


-3 


4.9 • 10' 


-3 


5.9 • 10' 


-3 


37 


8 


792 


94 


2310 


105 


445 


65 


8540 


282 


CH4 


8.3 • 10' 


-2 


9.7-10' 


-2 


1.1 • 10' 


-1 


32 


9 


346 


61 


570 


63 


379 


84 


1370 


157 


Ar 


4.3 • 10' 


-6 


4.9 • 10' 


-6 


5.8 • 10' 


-6 


31 


7 


333 


53 


557 


56 


288 


54 


1270 


104 


CO 


3.5 • 10' 


-7 


4.2 • 10' 


-7 


5.1 • 10' 


-7 


37 


8 


734 


88 


244 


102 


386 


59 


5840 


223 


C2H6 


7.5 • 10' 


-1 


7.6 • 10' 


-1 


7.7-10' 


-1 


3 


1 


50 


8 


96 


9 


17 


4 


108 


10 


CsHg 


7.3 • 10" 


-2 


7.4 • 10' 


-2 


7.5 • 10" 


-2 


3 


1 


50 


8 


96 


9 


17 


4 


108 


10 


C4H8 


1.4 • 10" 


-2 


1.4 • 10' 


-2 


1.4 • 10" 


-2 


3 


1 


50 


8 


96 


9 


17 


4 


108 


10 


HCN 


2.0 • 10" 


-2 


2.2 • 10' 


-2 


2.3 • 10" 


-2 


16 


4 


156 


50 


167 


56 


52 


13 


150 


52 


C4H10 


1.2 • 10' 


-2 


1.2 • 10' 


-2 


1.2 • 10' 


-2 


3 


1 


50 


8 


96 


9 


17 


4 


108 


10 


C2H2 


1.1 • 10' 


-2 


1.1 • 10' 


-2 


1.2 • 10' 


-2 


3 


1 


50 


8 


96 


9 


17 


4 


108 


10 


CH3CN 


9.7-10' 


-4 


9.9 • 10' 


-4 


1.0 • 10' 


-3 


3 


1 


52 


8 


112 


11 


17 


4 


109 


14 


CO2 


2.9 • 10" 


-4 


2.9 • 10' 


-4 


3.0 • 10' 


-4 


3 


1 


50 


8 


96 


9 


17 


4 


108 


10 




2.2 • 10" 


-4 


2.3 • 10" 


-4 


2.3 • 10" 


-4 


3 


1 


50 


8 


96 


9 


17 


4 


108 


10 



Tabic 3: Mole fractions of lakes species sensitivity to ground total pressure. 



Compound Molar fraction Molar fraction Molar fraction 

in liquid at in liquid at in liquid at 





Pjj^y = 1467 hPa 


P — Puuy 


X 0.9 




P — Puuy 


X 


1.1 




N2 


4.90 X 10-3 


4.07 X 10- 


-3 




-17% 


5.92 X 10' 






+21% 


CH4 


9.69 X 10-2 


8.26 X 10' 


-2 




-15% 


1.13 X 10' 


-1 




+17% 


Ar 


5.01 X 10-6 


4.29 X 10' 


-6 




-14% 


5.83 X 10' 


-6 




+16% 


CO 


4.21 X 10-^ 


3.49 X 10' 


-7 




-17% 


5.09 X 10' 


-7 




+21% 


C2H6 


7.64 X 10-1 


7.76 X 10' 


-1 




+1.6% 


7.50 X 10' 


-1 




-1.8% 


CsHg 


7.42 X 10-2 


7.53 X 10' 


-2 




+1.5% 


7.28 X 10' 


-2 




-1.9% 


C4H8 


1.39 X 10-2 


1.41 X 10' 


-2 




+1.4% 


1.37 X 10' 


-2 




-1.4% 


HCN 


2.09 X 10-2 (s) 


2.27 X 10' 


-2 


(s) 


+8.6% 


1.91 X 10' 


-2 


(s) 


-8.6% 


C4H10 


1.21 X 10-2 (ns) 


1.23 X 10" 


-2 


(ns) 


+1.7% 


1.19 X 10' 


-2 


(ns) 


-1.7% 


C2H2 


1.15 X 10-2 (ns) 


1.16 X 10' 


-2 


(ns) 


+0.9% 


1.13 X 10" 


-2 


(ns) 


-1.7% 


CH3CN 


9.89 X 10-^ (ns) 


1.00 X 10' 


~3 


(ns) 


+1.1% 


9.71 X 10" 


-4 


(ns) 


-1.8% 


CO2 


2.92 X 10-^ (ns) 


2.97 X 10' 


-4 


(ns) 


+1.7% 


2.87 X 10" 


-4 


(ns) 


-1.7% 


CeHg 


2.25 X 10-^ (ns) 


2.28 X 10' 


-4 


(ns) 


+1.3% 


2.21 X 10" 


-4 


(ns) 


-1.8% 
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Figure 1: Histogram of mole fractions of Titan's hydrocarbon lakes. A^* is the normalized 
number of mole fractions owning the value Xi computed within a given range of uncer- 
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Figure 2: Same as in Figure [T] 
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Maximum error considered in Monte-Carlo simulation (in fraction) 



Figure 3: Relative standard deviation a* corresponding to the mole fractions of N2, CH4, 
Ar and CO. a* arc plotted as functions of the maximum errors obtained for Pvap,i, ym,ii 
SHyap,i in Monte-Carlo simulations. Z^-'s have been set to 0.1 in all these computations. 
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Figure 4: Solid lines: time averaged ratios Ti/rcjHe from LEB08 models represented as 
a function of Titan's l atitude. Dash e d line s: the same ratios computed with the LAV09 
model and taken from lVuitton et al.l (|2008 ) in the case of CgHg. 



27 



